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ABSTRACT 

In the popular internal shock model for the prompt emission of gamma-ray 
bursts (GRBs), collisions between a series of relativistic shells generate lots of 
paired forward and reverse shocks. We show that the synchrotron emission pro- 
duced by the forward and reverse shocks respectively could peak at two quite 
different energy bands if the Lorentz factors of these two types of shocks are 
significantly different with each other (e.g., one shock is relativistic and the other 
is Newtonian). We then investigate whether this scenario is applicable to the 
case of GRB 080319B and find that a bimodal distribution of the shell Lorentz 
factors, peaking at ~ 400 and ~ 10 5 , is required. In addition, this scenario pre- 
dicts an accompanying inverse- Compton (IC) GeV emission with a luminosity 
comparable to (not much higher than) that of the synchrotron MeV emission, 
which can be tested with future Fermi observations. 

Subject headings: gamma rays: bursts — radiation mechanisms: nonthermal 



1. Introduction 



Since the pioneering works by Rees & Meszaros (1994) and Paczyhski & Xu (1994), it 
has been widely argued that the prompt emission of GRBs arises from internal shocks in a 
relativistic fireball that consists of a series of shells with different Lorentz factors and that 
the observed 7-ray emission is usually attributed to synchrotron or inverse- Compton (IC) 
emission from power-law electrons in the shocks (Meszaros & Rees 1997). Especially, in the 
case of IC 7-ray emission, we can also expect a bright synchrotron low-energy (e.g., optical) 
emission as seed photons for IC scattering. Based on this synchrotron and synchrotron self- 
Compton (SSC) emission scenario, Meszaros & Rees (1999) accounted for the bright optical 
flash of GRB 990123. 
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In this paper we put forward an alternative model, in which prompt optical and 7-ray 
emissions may be produced only by synchrotron radiation in internal shocks. It is natural 
to expect that internal shocks should consist of paired forward and reverse shocks, which 
are produced simultaneously by collisions between the shells in the fireball. In the case 
where the Lorentz factors of these two types of shocks are quite different, two groups of 
electrons accelerated by a forward and a reverse shock respectively during one collision 
are expected to reach different characteristic energies. Therefore, by combining the two 
synchrotron components contributed by these two shocks, we can get a bimodal photon 
spectrum, which might be able to account for some GRBs that have two spectral components. 
In particular, under some peculiar conditions, an optical peak can be found to fit the prompt 
optical emission of some GRBs such as the naked-eye GRB 0803 19B. 

For GRB 080319B, with its prompt 7-ray fluence (6.13 ±0.13) x 10~ 4 erg cm~ 2 (20 keV 
- 7 MeV; Racusin et al. 2008b) and redshift z = 0.937 (Vreeswijk et al. 2008), its isotropic 
equivalent 7-ray energy release is estimated to be E lyiso = 1.4 x 10 54 ergs (20keV — 7 MeV), 
which is among the highest ever measured. More surprisingly, this GRB was found to be 
associated by an extraordinarily bright optical flash peaked at a visual magnitude of 5.3 
(Racusin et al. 2008b), which is visible even for the unaided eye. Compared with the 
extrapolation of the 7-ray spectrum to the optical band, the observed flux density of the 
optical flash (~ 20 Jy) is about ten thousands times higher. In addition, thanks to the high 
time resolution, a fluctuating structure can be seen clearly in the optical light curve. 

For this rare mult i- wavelength prompt emission, any scenario invoking a single syn- 
chrotron component from only one emitting region is obviously not viable. Moreover, dif- 
ferent from some other bright optical flashes such as those of GRBs 990123 and 061007, an 
external shock origin for the optical flash of GRB 0803 19B is disfavored due to the short 
duration of the optical pulse and the lack of an increasing optical pulse duration throughout 
the whole prompt phase, as argued by Kumar & Panaitescu (2008). Instead, an internal 
dissipation origin may be plausible. As argued by Meszaros & Rees (1999) for GRB 990123, 
Kumar & Panaitescu (2008) and Racusin et al. (2008b) explained the prompt emission of 
GRB 0803 19B in the synchrotron and SSC emission scenario. However, as a direct con- 
sequence, they predicted a remarkably strong GeV emission, whose luminosity is about 10 
times higher than the observed MeV one. In contrast, for several bright GRBs detected by 
EGRET on CGRO, the GeV fluence is not higher than that in the MeV BATSE energy band 
(e.g. Sommer et al. 1994; Hurley et al. 1994). 

However, in our internal shock-produced two-component synchrotron emission scenario, 
the problem of GeV emission excess can be avoided, whereas the prompt optical and MeV 7- 
ray fluxes of GRB 0803 19B can be still explained if an unusually high variability of Lorentz 
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factors exists in the fireball. Moreover, since internal shocks can take place many times, 
we can naturally understand the fluctuating structure of the optical and 7-ray light curves. 
However, the emissions in the two energy bands are not strictly correlated with each other, 
since the durations of the shocks determined by the widths of the shells could be different 
and the specific shape of the light curves is dependent on the specific structure of the shells. 
For simplicity, these complications will be not considered in our simple model. 

Compared with the internal shock model suggested by Meszaros & Rees (1999) for GRB 
990123, our work differs from theirs in two aspects. First, as discussed above, we suggest 
that both optical and 7-ray emissions arise from synchrotron emission, differing from the 
synchrotron plus SSC scenario proposed by them. Second, in their works, only nonrelativistic 
to semirelativistic internal shocks are considered, which are produced by collisions between 
shells with similar Lorentz factors, while in our paper some relativistic internal (reverse) 
shocks would be involved in order to fit the observations of GRB 080319B. 

This paper is organized as follows: in §2, we describe the dynamics and synchrotron 
emission of internal forward-reverse shocks and the model parameters are expressed as func- 
tions of some observational quantities. In §3, we constrain the model parameters by the 
observations of GRB 080319B and then some implications from these results are discussed. 
In addition, the contribution to the prompt emission by IC scattering of the electrons is also 
considered with the inferred model parameters. Finally, conclusions and discussion are given 
in §4. 

2. Two-component synchrotron emission from internal shocks 

2.1. The dynamics and electron energy distributions 

In the internal shock model, the central engine of GRBs is assumed to eject a fireball 
consisting of a series of shells with different Lorentz factors 7 s heii during the prompt phase. 
Considering two shells that are ejected subsequently, for example, if the posterior shell 
(denoted by 4) moves more rapidly than the prior one (denoted by 1), a collision takes place 
at radius R is w 2^fc5t(f)~ 1 (Yu & Dai 2008) to produce an emission pulse, where St is the 
observed variability time and (fi z = 1 + z is introduced due to the cosmological dilution of 
time. Because of the collision, a pair of shocks (i.e., internal shocks) could arise: a forward 
shock propagating into shell 1 and a reverse shock propagating into shell 4. The shocked 
regions in shells 1 and 4 are denoted by 2 and 3, respectively. 

According to the jump conditions between the two sides of a shock (Blandford & McKee 
1976), we can calculate the comoving internal energy densities of the two shocked regions by 
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e2 = (721 - l)(472i + 3)nim p c 2 and e 3 = (734 - 1)(47 34 + 3)n 4 m p c 2 , where 721 or 734 is the 
Lorentz factor of region 2 or 3 relative to the unshocked region 1 or 4. The comoving proton 
number density rij of unshocked region i can be calculated by = L k ^/ (47ri?^7 J 2 m p c 3 ) for an 
isotropic kinetic-energy luminosity L k ^ and a bulk Lorentz factor 7^ (^> 1) of the unshocked 
shell. The mechanical equilibrium between the two shocked regions requires e 2 = e^, which 
yields 

(721 - l)(472i + 3) = n 4 = fjA 2 (1 \ 
(7 34 -l)(4734 + 3) m V74/ ' 

where L k> i = = L k is supposed. By assuming 74 ^> 71, the above equation leads to (Yu 
& Dai 2008) 

721 - 1 = 77 — + — - 1 « £ < 1, 734 = 77 — + — ~ 7T' 2 
2 V7 7i/ 2 V 74 7 / 2 7 

where 7 = 71(1 + y/2£) is the Lorentz factor of the shocked regions. This indicates that 
the reverse shock is relativistic and the forward shock is Newtonian. Moreover, following 
Dai & Lu (2002), the total number of the electrons swept-up by the forward and reverse 
shocks during a period of St can be expressed by N 2 = 2^/2^L k 5t/ (0 z 7im p c 2 ) and N 3 = 
L k 5t/ (0 z 7 4 m p c 2 ), respectively. 

Both forward and reverse shocks can accelerate particles to high energies and amplify 
magnetic fields. As usual, we assume that the energies of the hot electrons and magnetic 
fields are fractions e e and e# of the total internal energy, respectively. Thus, the strength of 
the magnetic fields is given 



,1/2 1 ( ^>z e B^k 



For the shock-accelerated electrons, a power-law energy distribution, due/d'ye oc 7~ p for 
7e > 7e,m> is assumed. The characteristic random Lorentz factors of these hot electrons in 
regions 2 and 3 are determined respectively by 

7e,m.2 = e e — (721 - 1) W h—i, (4) 
7e,m,3 = e e — 734 ~ 1 ~ ~ , (5) 

m e 2 m e 7 



1 In these expressions, the possible spreading of the shells (i.e. decreasing of L k ) is not taken into account, 
since a detailed description for the dynamic evolution is not necessary for our calculations. In this paper, 
the constant L k can be seen as an effective kinetic-energy luminosity. 
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where e e = e e (p — 2)/(p — 1). That 7 e ,m,3 ^ 7e,m,2, due to 74/7 ^> 2£, indicates that 
the characteristic energy of the reverse- shocked electrons is much higher than that of the 
forward- shocked electrons. Therefore, the resulting synchrotron photons emitted by these 
two types of electrons are expected to peak at two different energy bands and thus two 
distinct spectral components would be observed as in GRB 080319B. To be specific, the 
reverse shock is responsible for emission in a relatively high-energy band, while the forward 
shock contributes to a relatively low-energy component. 

In both shocked regions, the hot electrons with energies above 7 ei c,iW e c 2 lose most of 
their energies during a cooling time t C) j, where the cooling Lorentz factor is determined 
by 7e, c ,« — 67im e c(p z / (y i aTBf'jt c i ). The parameter y h defined as the ratio of the total 
luminosity to the synchrotron one, is introduced by considering the cooling effect due to the 
IC emission besides the synchrotron cooling. As pointed out by Ghisellini et al. (2000), 
the theoretical synchrotron spectrum arising from these electrons, calculated by using the 
standard assumption that the magnetic field maintains a steady value throughout the shocked 
region, leads to a spectral slope F u oc z/ -1 / 2 below ~100 keV, which is in contradiction to 
the much harder spectra observed. In order to overcome this problem, Pe'er & Zhang 
(2006) suggested that the magnetic field created by a shock could decay on a length scale 
(Ab,») much shorter than the comoving width (Aj) of the shocked region, i.e., \b,i = Aj/ fs,i 
(fB,i 3> 1). In other words, the shocked region can be roughly divided into a magnetized part 
immediately after the shock front and a further unmagnetized part. Under this assumption, 
the cooling time of the electrons should be determined by the time during which the electrons 
traverse the magnetized region, i.e., t c>i = St/ ' fs,i- Although the size of the magnetized region 
is reduced significantly by the field-decay effect (as found in §3.3), this region could be still 
wide enough for electrons to lose a great part of their energy when they traverse it. In this 
case, the cooling Lorentz factor 7 e c j of electrons is not much higher than 7 e , m ,j, so that the 
radiation efficiency of the electrons is not reduced drastically compared to the case without 
any magnetic field decay. 



2.2. Two-component synchrotron emission 

With the electron distributions and the magnetic fields described above, we can give 
the resultant synchrotron spectra using the method developed by Sari et al. (1998). The 
reference peak energies of the synchrotron spectra of the forward and reverse shocks are 
taken to be in optical and 7-ray bands, respectively. Then, the model parameters can be 
expressed as functions of some observational quantities. 

For the electrons in the magnetized reverse-shocked region, two break frequencies of the 
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synchrotron spectrum are given by 



<le 2 R 



2-Km e c(j) z 



(6) 



Qe 2 o 

^,3 = : 7e,c,3-^37 



2irm e cq 



4 V#wS44 

The peak flux density of the spectrum at z/ p = min[z/ m , z/ c ] reads 

<ft z m e c 2 a T N- 
^ 3 Audi 



(7) 




24irm p q e c 



where the parameter /# )3 is introduced because, at any moment, only a fraction 1//b,3 of 
the total reverse- shocked electrons locate at the magnetized region and other electrons in 
the unmagnetized region do not contribute to the synchrotron emission. The quantities in 
the left sides of equations (jEJ)-(|EJ) can be inferred from an observed prompt 7-ray spectrum, 
while the right sides are functions of the model parameters. We can therefore solve these 
equations to find the values of some model parameters, 



2.5 x 10 M erg s y 3 ,oe e ,ii 

X f^i,28^,P,3 ,-25 z/ c,3 2 20 I 

LI (9) 



^,3,20 



a 1 n4 2 -1/4 -1/4 .-3/4 
74 = 4 X 10 4 72 5^/3,0 £B,-l € eA 

X I OtO «L,28 ^y,p,3 -25^,3,20^,3,20 J > l 1U J 
^ R „m3 -4 7/4 3/4 .-3/4 

— 1/2 i r3/2 r 3/4 7/8 3/8 \ 

X I OTo ^,0.3^,28-^^,3,-25^,3,20^,3,20 J > l 1J -J 

where and hereafter the convention Q = lCPQa, is adopted in cgs units. The quantities in 
the brackets are basically determined by the observational data and the values of L^, 74, 
and /b 3 are modulated by the remaining free parameters. In particular, 74 and /b 3 are 
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strongly dependent on 7 that can be constrained by optical observations. The fact that 
is independent of the parameter f B 3 indicates that the hypothesis of magnetic field decay 
does not increase the energy requirement of the model. This is because all the reverse- 
shocked electrons, when they traverse the tiny magnetized region at different times, have 
released a great part of their energy to 7-rays (as indicated by f m ,3 ~ ^,3) via synchrotron 
emission. The observed 7-ray emission is mainly contributed by the emission from this tiny 
region behind the shock front, no matter whether the other part of the shocked region is 
magnetized or not. 

In order to study the properties of the forward shock, we now focus on a possible low- 
energy emission, whose spectral information is however not as rich as the 7-ray component. 
We calculate the peak frequency of the synchrotron spectrum produced by the forward shock 

as 



In order to explain the optical spectrum measured by Raptor after 80 s for GRB 0803 19B 
(Kumar & Panaitascu 2008; Wozniak et al. 2008) in the next section, we assume that 
this peak frequency is below the optical band (i.e., v m p < v Q = 5 x 10 14 Hz), which yields 
7 > 7 = 288 yl^e^^e^tiStQ 1 ^ 2 L*^^. This requirement is however not necessary for common 
GRBs. For z/ m2 < v < z/ c 2 , we can calculate the optical flux density and the synchrotron 
self-absorption thickness at u Q by 



Qe 



2imi e c(f) z 




(12) 




(13) 




(14) 



where 



F, 



u,p,2 



(f) z m e c 2 a T N 2 
And 2 L 3q e f B ,2 



V7rn e a T ( 4> 2 z t B L\e 



) 



1/2 



6%m p q e c 3 / 2 y rf|/l,27' 



(15) 



5q e N 2 



r. 



sa,m 



5 27e 5 , m , 2 47^^,2 
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5y/7m 5 e q e 

28^5/2 I e B efp B ^ 



1/2 



4tl0 



Combining equations (|T2l) . ({TBI and ffT4"l) . we write the parameters 7 and /b,2 as 
of 6b and e e , 



(16) 
functions 



7 



420 u 1/16 e 1/16 r 1/16 

v ( Ay- -5 / 8 //)^ 1 / 4 ^ 1 / 2 P 1 / 4 r*l/16 -1/4 
I ot Tz,0.3 a L,28 r u,o,-22 lj k,52A 1 sa,o,-l 



(17) 



B,2 



180 y 3 



3(6+p)/16 (2+3p)/16~-(34-13p)/16 



Z B,-1 



-e,-l 



c,(14+p)/8 ,(6+p)/4 .-(6+p)/2 
ot ^2,0.3 U L,28 



v p -(6+p)/4 r *3(6+p)/16 (2+ P )/4 
A r i/,o,-22 -^fc.52.4 'sa,o,-l 



where the kinetic-energy luminosity Lk = 2.5 x 10 52 erg s _1 (?/3/e e )L^ 52 4 obtained in equation 
([9]) has been substituted. The value of 7 is mainly determined by the observational data 
and insensitive to the remaining free parameters e# and e e . For common GRBs, the value 
of T sa>0 is thought to be very high (e.g., ~ 10 3 ) and thus the value of 7 could be lower than 
100, which leads to hundreds for 74. However, for GRB 080319B, r S(l o is deemed to be not 
larger than unity to ensure a bright optical flash and we suggest r saj0 = 0.1 as a reference 
value hereafter. 



3. Application to GRB 080319B 

GRB 080319B triggered the Swift-Burst Alert Telescope (15-350 keV) at T =06: 12:49 
UT on March 19, 2008 (Racusin et al. 2008a) and was simultaneously detected by the Konus 
7-ray detector onboard the Wind satellite (20 keV-15 MeV; Golenetskii et al. 2008). The 
time-averaged Konus- Wind 7-ray spectrum can be fitted well by the Band function (Band 
et al. 1993) with a low energy slope of 0.855^^3 below the peak of E p = 675 ± 22 keV and 
a high-energy slope of — 3.59^g'g2 above the peak (Racusin et al. 2008b). The burst had a 
peak flux of F p = (2.26 ± 0.21) x 10~ 5 erg cm~ 2 s _1 and thus the peak isotropic equivalent 
luminosity was L Pjiso = (1.01 ±0.09) x 10 53 erg s -1 . Using the values of F p and E p , we roughly 
estimate the peak flux density of the 7-ray spectrum, F p j E p w 14 mJy Compared with the 
extrapolation of the 7-ray spectrum to the optical band, the observed flux density of the 
optical flash (~ 20 Jy) is about ten thousands times higher. 
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3.1. The model parameters 

Adopting z = 0.937 (d L = 1.9 x 10 28 cm), 5t ~ 3 s, F„ )0 ~ 20Jy (at z/ D = 5 x 10 14 Hz), 
Fv,v,z ~ 14mJy, fof mi 3 ~ 675keV and denoting z/ Cj3 = a; x z/ mi3 = 10°xo x u m ^ for GRB 
0803 19B, we derive the model parameters, 

L k ~ 2 x 10 54 erg s- 1 ^ /2 y 3 ,od, (19) 
7 ^ 400 xrysT^iC?". ( 20 ) 
74 - 9 x 10 4 x?l™vty%\!XCl ( 21 ) 

/ B>2 ~ 700 x f + P )/3^3(6 +P )/16 e g + ^)/16.-(34-13 P )/16 ; (22) 

/ fl>3 - 7 x io 3 4 /4 y 3 3/ o 4^C?> ' (23) 

where a fiducial value of unity is assumed for y 3 , which will be proved in §3.2. Besides 
a constraint by the maximum allowed equipartition value (e e < 0.3 and e# < 0.3), the 
remaining free parameters e# and e e satisfy e e < 0.09(xy|e^) -3 / 26 given 7 > 7. The upper 
limit of e e is insensitive to the value of e# (strictly, with a decrease of €3, the upper limit of 
e e increases slightly). Taking e e < 0.09 (e e < 0.3) as a conservative estimate, we find: 

(i) L fe > 2 x 10 54 erg s _1 . This is a natural result due to the high observed 7-ray luminos- 
ity (~ 10 53 erg s _1 ) of GRB 080319B. The MeV 7-ray radiation efficiency of the reverse shock 
can be estimated by rj ~ 0.05x 1 ^y^ ) \e e -\. We next calculate the total isotropic-equivalent 
energy release of GRB 0803 19B, 

E Kiso = 2E lMo /r] ~ 5 x 10 55 x 1/2 y 3 ^-i_ i ergs, (24) 

where a factor 2 is introduced by considering a similar amount of energy carried by the 
forward shocks. Using a very small jet angle 9j = 0.2° that is found by Racusin et al. 
(2008b), we get the beaming-corrected energy release of GRB 080319B, 

E kJet = E k , so 6]/2 ~ 10 51 x 1/2 y 3 ,od ergs, (25) 

which is a typical value for common GRBs. 

(ii) 7 ~ 400. This is a typical value for the Lorentz factor of merged GRB ejecta after 
internal shocks. Thus, we can obtain the internal shock radius for GRB 080319B, 

R ls » 2 7l 2 c(5t0; 1 » 10 16 cm fto-s^aaT^e, (26) 

which is larger than those of common GRBs (~ 10 14 cm). So, the dynamic influence of the 
circum-stellar medium on the leading shell (71 ~ 7 ~ 400) needs to be assessed. We calculate 
the deceleration radius of the leading shell via 

Rdec = ^2l c 2 ~ 5 X 10l6cm ^-1^0.5^0.37i:2 2 6^,54.3, (27) 
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where a wind-like circumburst medium (p = Ar~ 2 with A = 5 x 10 11 g cm _1 A*) is assumed. 
Moreover, as claimed by Racusin et al. (2008b), a tenuous wind with an upper limit of 
A* ~ 0.03 is required by the afterglow data of GRB 080319B. In this case, the deceleration 
radius is larger. Therefore, we conclude that the deceleration of the shells can be ignored at 
the times where internal shocks among the shells occur. 

(iii) 74 > 10 5 . This high Lorentz factor is allowable for acceleration of an initial fireball 
with very low baryon contamination (Piran 1999). 

(iv) fs,3 > 10/b,2- Our constraint on fs,3, which is much larger than unity and much 
less than A/A s (\ s is the plasma skin depth), is consistent with that found by Pe'er & Zhang 
(2006) for other GRBs. However, the physical underpinning of these values of J'b is unknown. 
So the difference in fs for different shocked regions lacks a reasonable physical explanation 
and it might be related to different shock strengths of the forward and reverse shocks. 

From the above results, we can see that most of the inferred model parameters are 
reasonable and acceptable even in common GRBs, except for an unusually high variability 
of Lorentz factors denoted by 74/71 ~ 300. Although the 7-ray emission of GRB 080319B 
seems to be not unusual, the relatively high value of St ~ 3s (versus ~10 ms for com- 
mon GRBs) implies an unusually large internal shock radius, which ensures the synchrotron 
self-absorption frequency below the optical band and reduces the magnetic field strength. 
Therefore, in order to produce sufficiently energetic 7-ray photons, it is necessary to invoke 
some highly relativistic internal shocks that require high variability of Lorentz factors in the 
fireball. Since the GRB central engine is far from being thoroughly understood, it is difficult 
to demonstrate whether the central engine can produce such a drastically varying fireball or 
not, but some possible origins can be still imagined. In the collapsar model, for example, a 
shell passing through the envelope of a progenitor star should sweep up and clear away the 
envelope material, leaving a channel behind the shell. A following shell will pass through 
this clear channel and thus have a very low baryon contamination. This might lead to a 
highly relativistic (7 s heii ~ 10 5 ) shell. However, due to lateral diffusion of the channel wall, 
this channel will be possibly contaminated by baryons again some time (~ order of seconds) 
later. Subsequently, such a switching-on-and-off process of the channel repeats again and 
again. As a result, relatively slow and rapid shells are generated alternately. In reality, this 
process is unlikely to be so regular because a slow/rapid shell could be followed by another 
slow/rapid shell. Therefore, the temporally-correlated optical and 7-ray emissions from this 
process could be polluted by the emission due to the collisions between slow-slow or rapid- 
rapid shells. This along with other effects (e.g., different shock-crossing times of shells and 
so on, see §4) may lead to the fact that the prompt optical and 7-ray emissions are not 
correlated finely, as observed in GRB 080319B. 
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To summarize, if some unknown physical processes of the central engine can give rise to 
a fireball in which the Lorentz factors and densities vary drastically in a variability timescale 
of few seconds, our internal shock-produced two-component synchrotron emission scenario 
may account for some fundamental features of the optical and MeV 7-ray emissions of GRB 
080319B. 



3.2. Inverse Compton emission 

We denote r« = axNi/ (4:nR 2 fB,i) and t\ = utNi/ (AtiR 2 ) as the Thomson optical depth 
of the magnetized and unmagnetized regions, respectively. The SSC emission from the two 
magnetized regions is considered first. Following Sari & Esin (2001) forp = 2.5, we estimate 
the two break energies of the SSC spectrum contributed by the forward shock by 

hl/ ssc_ 2 hl . _ 1 ^ w J/ 16 ,, 3 / 8 , 3 / 8 ? 29 / 8 (oq\ 

h^f = H,c,2^c,2 = 70 GeV xly^ y yi_ v (29) 
The peak flux at v^P can ^ e estimated by 

r,,piSSC ,,SSO P /,,SSC /,,SSC\ "(P" 1 )/ 2 

~ 2 x 10" 8 erg cm^s" 1 (30) 

For the reverse shock, as the SSC peak enters the Klein-Nishina regime, the real peak energy 
can be determined by (Gupta & Zhang 2007; Fragile et al. 2004) 

»■«?, = ^ = 60 GeV (31) 
at which a negligible flux is found from 



,.ssc 



1/3 



27^,3^,3 1 

~ 8 x 10~ n erg cm^s -1 . (32) 



According to the definition of parameter y iy we obtain 



rSSC \yj? 1SSC 

Ite = l + ^<l + fSr-«l-2, (33) 



opt 



rSSC \ v p 1SSC 

^3 = 1 + 7^1 + ^^1, (34) 
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where [^i^] opt = 10 _7 erg cm~ 2 s~ 1 and [vF v }mcV = 2.3 x 1CT 5 erg cm -2 s -1 . This indicates 
that the SSC emission is relatively unimportant and the estimations in equations ffT9l) - (l32l 
under the assumption, y 2 ~ 1/3 ~ 1, are self-consistent. 

Although a part of the electron energy can be released via synchrotron and SSC emis- 
sions when the electrons traverse the magnetized regions, about half energy of the reverse- 
shocked electrons and almost all energy of the forward-shocked electrons are still held by 
the electrons when they enter into the unmagnetized regions. This remaining energy can be 
no longer released via synchrotron radiation because of the lack of magnetic fields, but can 
via external inverse Compton (EIC) scattering due to the existence of the radiation fields. 
Determined by this EIC cooling, the cooling Lorentz factor of the electrons in both of the 
unmagnetized regions reads 



m P c(i>? 



where the radiation energy density contributed by the synchrotron radiation from the two 
magnetized regions can be calculated by w 7 = (y 2 + 2/3 — 2)(B 2 /8ti) = Y\B 2 /8ti) by consid- 
ering B2 = -B3 = B. Then we have 

f _ Q7rm e C(j) z _ H/32 yt -l -11/16 -H/16.11/16 ( ^ 

7e ' c ~ Ft aT 5 2 7 5t ~ r -o-7^3,o e B ,-i ^-x ■ W 

and 7g C < (7e,m,2, 7e,m,3)- F° r an EIC spectrum produced by upscattering the seed photons 
from magnetized region j by the electrons in unmagnetized region i, we calculate its two 
characteristic break frequencies by 



v 



U) 

L.i 



H,L^L d , (37) 



41 = Kh^hj, (38) 
and estimate its peak flux at the peak frequency z/j^ roughly by 

K.lS - ^.r}F u>P!(j) (v^/vty HP ~ 1] ' 2 , (39) 

where the subscription L represents the low break frequency of the seed photons and the low 
break Lorentz factor of the target electrons, while H represents the high ones. Considering 
the two synchrotron components for seed photons and the two population unmagnetized 
electrons, four EIC components are expected: 

(i) For i = 3 and j = 2, we have 

= 27£>m,2 = 0.5 keV *o 1/2 ^o 2 7 %>S^1,-i, (40) 
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h^H )3 = 27^,3^^,2 = 6 TeV xf /32 y 2 ^y% 16 e B ^^if. (41) 

(2) 

Although the peak frequency v H 3 is just around the Klein-Nishina break determined by 
hu^ N3 = 7 2 mgC 4 /(/iz/ c 2 ) ~ 6 TeV, we may still derive a relatively low flux from equation 

i vF »]pl ~ 8 x 10_8 er S cm-V 1 . (42) 

(ii) For i = 3 and j = 3, we have 

= 2 7e > c , 3 = 0.3 GeV xl^ 7 y^%%^% (43) 

M? 3 = 27 e 2 m ,3^^,3 = 0.6 EeV x^y^e^^. (44) 

(3) (3) 

The peak frequency u H 3 is however higher than the Klein-Nishina break frequency, hv KNZ = 
r y 2 m 2 e c i I \hp m ^) ~ 60 GeV, and thus the real peak flux can be calculated at v^nz ^° ^ e 

[vF v ]®l ~2x 10~ 7 erg cnrV 1 . (45) 

Summarizing the above two cases, the EIC process of the reverse-shocked electrons produces 
a TeV and a GeV emission component, both of which are much weaker than the observed 
synchrotron MeV emission with a flux of ~ 2.3 x 10~ 5 erg cm~ 2 s _1 . 

(iii) For i = 2 and j = 2, we have 

K 2 2 = 27j> mJ = 0.5 keV x^YI^y^-A^ (46) 

**& = = 10 WV xf^lvT4^C". (47) 



(2) 

and the peak flux at i/# 2 



~ 7 x 10_7 er S cm _2 s _1 . (48) 



Due to the stronger synchrotron MeV emission, this relatively weaker MeV emission com- 
ponent is likely to be covered up. 

(iv) For 2 = 2 and j = 3, we have 

= 2 T t> w>2 = 0.3 GeV 4 /1 Vr 2 7 2/ 3 :o 11/8 ^-( 8 < 1 5, (49) 

Mgj = H,m,Mm >3 = 0.9 GeV (50) 
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(3) 

Distinct from the emission components discussed in the previous three cases, the flux at v H 2 , 

[vF v ] { *l ~2x 1(T 5 erg cnrV 1 , (51) 

indicates a strong GeV emission that is as strong as the synchrotron MeV emission. This 
means that the energy of the forward-shocked electrons would be mainly released by upscat- 
tering the synchrotron MeV 7-ray photons from the reverse shock, while the reverse-shocked 
electrons lose a great part of their energy via synchrotron cooling directly. Due to the strong 
emission within the energy regime from sun-GeV to GeV, electron-positron pairs might be 
produced by collisions between the sub-GeV and GeV photons. According to the results 
obtained above, we can give approximately an upper limit luminosity of Li im ~ 10 53 erg s _1 
for the sub-GeV (e 7 ~ O.lGeV) emission. Then, the optical depth due to pair production 
interactions can be roughly estimated as 



T. 



™ ~ \ aT J^E^i ~ °- 3 Llim ' 537 ^r 6 \ (52) 



which indicates that the pair production effect is not significant. For simplicity, the further 
contribution from the secondary electrons is ignored here. 

To summarize, the contributions by the SSC and EIC emission to the observed optical 
and MeV 7-ray emissions are insignificant and the two synchrotron components are dom- 
inant in the observed bands for GRB 080319B. In contrast, some higher-energy emission 
components would be produced by the EIC process. Although most of these components are 
weak, the flux of the strong GeV emission can reach as high as ~ 10 _5 erg cm~ 2 s _1 , which 
is comparable to the observed MeV one. 



4. Summary and conclusions 

In the popular internal shock model for the prompt emission of GRBs, paired forward 
and reverse shocks are produced by collisions between some relativistic shells with different 
Lorentz factors. In this paper, we have considered this model in an unusual situation where 
a bimodal distribution of the shell Lorentz factors exists. As a result, the Lorentz factors 
of the forward and reverse shocks are quite different (i.e., the forward shock is Newtonian 
and the reverse shock is relativistic) and the resulting two-component synchrotron emission 
is expected to provide a new scenario for some seldom GRBs that have two spectral peaks 
in the prompt emission. 

As an example, we compare our scenario with the recently-observed naked-eye GRB 
080319B and constrain the model parameters by fitting the observations. We find that, 
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on one hand, the optical and MeV 7-ray fluxes of this unique GRB could be explained in 
our two-component synchrotron emission scenario, if some unknown physical processes of the 
central engine (e.g., the picture described in §2.1) can give rise to a fireball where the Lorentz 
factors of ~ 400 and ~ 10 5 appears alternately in a variability timescale of few seconds. On 
the other hand, although the internal shock-produced emission can roughly account for the 
fluctuating structure of the light curves and the mild temporal correlation between the optical 
and 7-ray emissions, it is still difficult to explain clearly why the observed optical emission 
varies relatively slower than the 7-ray emission. We speculate that this difference between 
the light curve variabilities may be due to a complicated distribution of Lorentz factors in 
the fireball, an inhomogeneous structure of each fireball shell, different shock-crossing times, 
and other more realistic properties of the system. So, a more detailed simulation is required 
to improve our present model. 

Finally, for high energy emission, the synchrotron plus SSC scenario suggested by Kumar 
& Panaitescu (2008) and Racusin et al. (2008b) predicts significant GeV 7-ray emission by 
considering the second order IC-scattering, the flux of which is about 10 times higher than 
the observed MeV one. In contrast, our model predicts a relatively weaker GeV component, 
whose flux is lower than or at most comparable to that of the synchrotron MeV emission. 
Therefore, future observations for high energy counterparts of GRBs by the Fermi Space 
Telescope are expected to be able to discriminate these two models. 
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